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Ejercicio 1. Sea g un campo vectorial, cuyo campo de rotores es h (x, y, z) )= (z-2x, y, x 2 + y 2 + z), y sea la superficie S = 
{(x, y, z)eM 3 : x 2 + y 2 + z 2 = 25, y > - 4). Calcular la circulación de g sobre la curva C, borde de S, indicando en un gráfico 
la orientación asignada a C. 



Ilustración 1. Curva C y superficie S correspondientes a los objetos geométricos del enunciado 

El gráfico anterior muestra los objetos geométricos del enunciado: la superficie S y su contorno C, con una cinta negra. Dado 
que g es 1 II C 00 (M 3 ), que C es un lazo simple regular, como lo prueba la parametrizaciórfii y: [0, 2 tt] — » M 3 tal que y(t) = (3 
cos(t), -4, 3 sen(t)), con velocidad y': [0, 2 ti] — » M 3 tal que y'(t) = (-3 sen(t), 0, 3 cos(t)), que orienta la curva como lo 
indica la figura siguiente. Ahora, si definimos la superficie^) = ={(x, y, z)eM 3 : x 2 + z 2 < 9, y = -4}, disco que puede 
parametrizarse en JZ = [0, 3] x [0, 2n) como So = T(JZ), con T(r, t) = (r cos(t), -4, r sen(t)), campo de normales dado por 
el producto vectorial Ñ(r, t) = T r (r, t) x T t (r, t) = (0, -r, O) 111 IV . Con las orientaciones así planteadas iv , e indicadas en la 
ilustración 2, podemos aplicar el teorema de Stokes para calcular lo pedido, cuyo valor es 36 ti. 




ñdS = 



h(T(r, t)) • N(r, t)drdt = 4 Jj r dr dt = 


2n 

= 4 I dt I r dr — 36n 


La anterior resolución es muy ineficiente, la parametrización del disco So bien pudo evitarse, dado que se trata de la 
porción de un plano con versor normal inmediato (constante en todos los puntos de la superficie) y está dado por: 
n— (0, — 1, 0), de modo que, para todo punto (x, -4, z) en So se tiene que /¡(x, — 4, z) • n(x, — 4, z) =4, de modo que 

JJ h-ñdS = 4 JJ dS = 4 área(S) = 4(97i) = 36n. La habilidad de escoger la mejor solución se adquiere desde la 
práctica, y comparando más de una resolución. 

Una resolución todavía peor que la primera consiste en ignorar el Teorema del rotor, y recuperar un campo vectorial g a 
partir de su campo de rotores (un problema cuya solución está determinada salvo el gradiente de un campo escalar v ): uno 


I Su campo de rotores es C 00 en R 3 , de modo que debe serlo él mismo. Se dice que g es un potencial vectorial del campo de rotores h = rot( g), lenguaje 
similar al utilizado cuando se dice que cp es un potencial escalar del campo vectorial de gradientes f= Vcp. También se dice que h es el rotacional de g. 

II Es claro que y es C 00 en su dominio (sus componentes son senos y cosenos), y además si y(t¡) = yiif), esto es, si (3 cos(ti), -4, 3 sen(ti)) = (3 cos(t 2 ), -4, 3 
sen(t 2 ))es necesariamente que cos( ti ) = cosOA, y sen( t¡ ) = sen(t 2 ) y de aquí es t¡ = t 2 a menos que difieran en 2n, lo que sólo sucede en los extremos del lazo. 
Además, para ningún valor de t el vector y’(t) = (-3 sen(t), 0, 3 cos(t)) es nulo, pues los valores que anulan la función seno no anulan la función coseno. 

III Es inmediato que esta parametrización es C°°CR), y que la normal sólo se anula en el subconjunto Hn ={(r, t) e R : r = 0}, que es de contenido nulo en R , de 
modo que se trata de una parametrización regular en casi todas partes en R. Por otra parte, la parametrización es inyectiva en casi todas partes (creo un 
ejercicio muy instructivo la prueba de esta afirmación, siempre que lo haga uno mismo). 

IV Las orientaciones se hallan concordadas como lo exige el teorema de Stokes, que exige preservar la orientación del espacio introducida por la terna de referencia 
en la que se calcula el rotor (en toda nuestra asignatura, mantenemos el sistema dextrógiro o inglés). 

v Puede verse en Santaló, L. (1993). Vectores y tensores con sus aplicaciones (Décimocuarta edición ed.). Buenos Aires: Eudeba, Capítulo V, Fórmula Integrales, 
§23.3, p.178, teorema 4, con su demostración constructiva de uno de tales campos. También, en Apóstol, T. (1980). Calculas volumen 2. Cálculo con funciones de 
varias variables y álgebra lineal, con aplicaciones a las ecuaciones diferenciales y a las probabilidades (Segunda edición en castellano [Original: Calculus II, 
Multi-variable calculus and linear algebra, with applications to differential equations and probabilty] ed., Vol. II). (F. Vélez Cantarell, Trad.) Barcelona: Reverté, 
Capítulo 12, Integrales de superficie, §12.6 Reconstrucción de un campo vectorial a partir de su rotacional, Teorema 12.5, p.547; es recomendable leer el ejemplo 
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de tales campos es g(x, y, z) = (0, x 3 /3 + y 2 + zx, -x y + zy), y todo otro campo g* es g* = g + Vcp para algún campo 
escalar cp. La integral curvilínea de g* es igual a la de g, puesto que la circulación del Vcp es nula sobre C (¡es un lazo!), de 

modo que podemos calcular de modo directo / g • di, con lo que añadiremos al cálculo ya hecho para recuperar g, el de 
esta integral curvilínea (hacerlo, como ejercicio, para concluir que sería una muy mala elección). 



Ilustración 2. Las orientaciones adoptadas para la superficie S G y la curva C que es su contorno 


Este comentario se origina en la observación de que, a menudo, el teorema del rotor es desperdiciado, y hasta 

convertido en perjudicial. La potencia del teorema proviene de que la circulación del campo sobre C iguala el flujo de su 
rotor sobre cualquier superficie suficientemente suave y orinada adecuadamente que tenga como contorno la curva C; es 
justamente esa la libertad que deja lugar al que enfrenta el problema, a seleccionar 'la mejor' en algún sentido que facilite 
las cosas. En la resolución anterior hemos utilizado el disco plano, pero no han faltado exámenes en que se ha intentado 
hacerlo con la superficie S, o con el casquete esférico en y < - 4, lo que es una forma asegurarse serios problemas. 
Exagerando, podríamos imaginar una persona eligiendo (legítimamente) la superficie del siguiente gráfico: lo legítimo no 
coincide con lo conveniente. 



Ilustración 3. Una buena manera de convertir a Stokes en perjudicial 


de un campo solenoidal que no es un rotacional, para remarcar que la condición necesaria y suficiente que enuncia el teorema se asegura en abiertos simplemente 
conexos. Finalmente, el mismo punto es tratado en Rey Pastor, J., Pi Calleja, P., & Trejo, C. (1968). Análisis Matemático II. Cálculo infinitesimal de varias 
variables. Aplicaciones (Séptima edición ed., Vol. II). Buenos Aires: Kapelusz., §91.6, c (otras relaciones), p.518. 
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Ejercicio 2. Sea heC^M * * 3 ) un campo escalar y feC^M 3 ) un campo vectorial solenoidal (div f idénticamente nula en su 
dominio), tal que para todo punto P = (x, y, z) de M 3 es Vh (P) _L f(P). Sea S = {(x, y, z)el 3 : x 2 + y 2 + z 2 = 4}. Probar que el 
flujo del campo vectorial g(P) = h(P) f(P) a través de S es nulo. 


El ejercicio resulta casi inmediato si se tiene presente la importante identidad vectorial vi V-(h f) = Vh • f + h V- f, que con 
las hipótesis del ejercicio indica que div ( g ) = div(h f) = V-(h f) = 0. Siendo la superficie S una esfera, es claro que es 
regular y orientable, y entonces puede orientarse con normal saliente vii . Por otra parte, el campo vectorial g, por definición, 
es C 1 (M 3 ), dado que es el producto de un campo escalar y uno vectorial de ese grado de suavidad, por lo tanto está definido 
en un abierto simplemente conexo (esto es M 3 ) que contiene a S y a JVC = {(x, y, z)ei 3 : x 2 + y 2 + z 2 < 4}. Aplicando 
entonces el teorema de la divergencia viii resulta: 




div(hf^dV — 0 


El cálculo anterior prueba que, efectivamente, el flujo del campo vectorial g a través de S es nulo, 
resultado subsiste reemplazando la esfera del enunciado por cualquier superficie S frontera de un macizo JVC, dado que el 
argumento esencial depende del carácter solenoidal del campo que fluye a través de esa superficie, y no de la forma de S. 
Las identidades vectoriales señaladas en el pie de página anterior (¡21 identidades!) permiten reducir la complejidad de un 
problema al mínimo. 


a 


Resulta claro que el 


Ejercicio 3. Sea f: M 3 — >M 3 tal que f(x, y, z) = (x, y/2, z), y C la curva parametrizada por la función vectorial y: M -> M 3 tal 
que y(t) = (e*, e 11 / 2 , e 11 ). (a) Comprobar que C es la línea de campo de f que pasa por Po = (1, 1, 1). (b) Hallar el flujo de f a 
través del plano normal a C en Po en el primer octante, orientado con el campo de normales alejándose del origen de 
coordenadas. 


(a) Siendo /(y(t)) = (e í ,-e í/2 ,e í ), mientras que y (t) = (e l ,-e t/2 , es tal que Po = y(0), resulta inmediato que vale la 

expresión /(y(0)) = /(P 0 ) = ( l) = y'(0), lo que prueb^b pedido ix . 

(b) El plano normal a C en Po es normal al vector N = (2, 1, 2) y X está en tal plano sii (X-Po)-N = 0, esto es que su ecuación 
es 2x + y + 2z = 5. Luego la superficie S es S = {(x, y, z)el 3 : 2x + y + 2z = 5, x>0, y>0, z>0}, y el cálculo del flujo puede 
hacerse mediante la parametrización T en el recinto JZ = {(u, v)eM 2 : 0<v<5, 0<u^ 1/2 (5 - v)}, siendo T: JZ ^M 3 tal que 

T(u, v) = (u, v, 5/2 - u -v/2 ), con N(u, v) = T u (u, v) x T v (u, v) = (1, % , 1), de modo que se tiene S = T (JZ), con la 
orientación pedida. De esta manera, resulta que 



f-ñdS = Jj f(T(u,v))-N(u,v)dudv = 



( v 5 

\ 2 , 2 ~ U ~ V/2 


1 \ 1 f 5 f—, N 625 

l,-,l]dudv = -\ dv I (10 — v)du = — — - 

Z / T" JQ Jq T'O 


V1 La identidad vale con heC/R 11 ), feC/R 11 ), para todo n natural (para n = 1 es la conocida expresión de la derivada de un producto, de la que esta identidad es una 
generalización). La prueba es inmediata, llamando a la derivada respecto de la i-ésima componente, ( f )¡ al i-ésimo campo escalar de f, se tiene la siguiente 

cadena de igualdades div(hf ) = X? =1 A(A/)¿ = E?=i [A(A)(/)¿ + h P (/)¿1 = 'Ta=iD i (Ji)(f) i +'Tn=ihD i (f') i = V(/t) •/+ hdiv(f), la primera igualdad por 
definición de divergencia, la segunda por aplicación de la derivada del producto de funciones escalares, la tercera por asociatividad de la suma de funciones 
escalares, la cuarta por definición de gradiente, producto escalar y divergencia. En el excelente y sensato texto Marsden, J. E., & Tromba, A. J. (1991). Cálculo 
Vectorial (Cuarta edición en español del original Vector Calculus, Third Edition. ed.). (M. López Mateos, & S. Adarve, Trads.) Wilmington, Delaware, EUA.: 
Addison-Wesley Iberoamericana, Capítulo 3, Funciones con valores vectoriales, §3.5. Cálculo vectorial diferencial, tabla 3.1 (algunas identidades comunes en el 
análisis vectorial), p. 231, se encuentra esta identidad (es la número 8) junto a otras veintiuna; la identidad 8 está probada para el especial caso de R 3 en el ejemplo 1 
de la página232. Finalmente, también puede verse en Rey Pastor, J., Pi Calleja, P., & Trejo, C. (1968). Análisis Matemático II. Cálculo infinitesimal de varias 
variables. Aplicaciones (Séptima edición ed., Vol. II). Buenos Aires: Kapelusz, §91.6.b, expresión [91.36]. 

vu Una parametrización C" en R = [0, ti] x [0, 2n) que logra S = T (R) es f (cp, 0) = (2 sen (cp) eos (0) , 2 sen(cp)sen ( 6 ), 2 eos (<p)), que orienta S y es regular 
en todo “R , excepto los conjuntos Hi = {(cp, 0)e^ : cp = 0}, FE = {(cp, Q)eR : cp = ti}, de contenido nulo en R y cuyas imágenes a través de T son, 
respectivamente, el polo norte y sur de la esfera. En tales puntos (y sólo en ellos) 6 ) x T g ( cp , 6) = 4 sen ( cp ) {sen ( cp ) cos(0) , sen(cp)sen ( 6 ), eos (<p)) es 

el vector nulo. Además, T es, exceptuando Hi y EL, inyectiva en R (la prueba de esta afirmación es instructiva). 

vm Admitido conocido el resultado de que, para un campo vectorial en R 3 , de clase C 1 en un abierto simplemente conexo que contiene a una superficie cerrada S 
regular por partes, simple y orientada con normal saliente, el flujo del campo a través de S coincide con la integral triple de su divergencia en la región M de la que 

5 es frontera. Las hipótesis pueden debilitarse, y la dimensión del espacio ampliarse, resultando teoremas más generales; puede verse una presentación introductoria 
en Lang, S. (1990). Introducción al Análisis Matemático (Primera edición. [Original 1968, Undergrauate Analysis] ed.). (M. López Mateos, & M. Muñoz Mella, 
Trads.) Wimington, Delaware.: Addison Wesley Iberoamericana, Capítulo 20, Formas diferenciales, pp. 442-458. 

1X Se acepta conocida la definición de que y(t) es una línea de flujo del campo vectorial /(P) en el punto P = y(t) sii /(y(t)) = que puede verse en, por 

ejemplo, Marsden, J. E., & Tromba, A. J. (1991). Cálculo Vectorial (Cuarta edición en español del original Vector Calculus, Third Edition. ed.). (M. López Mateos, 

6 S. Adarve, Trads.) Wilmington, Delaware, EUA.: Addison-Wesley Iberoamericana, Capítulo 3, Funciones con valores vectoriales, §3.3, Campos vectoriales, 

p.216, definición. 
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Seguramente, el modo más sencillo de resolver este ejercicio consista en comprobar a simple vista que el flujo a través 
de los tres triángulos definidos en cada plano coordenado por las trazas del plano y los ejes coordenados es nulo, y como la 
divergencia del campo es constante (e igual a 5/2), el flujo es (por el teorema de la divergencia) esa constante por el 
volumen del tetraedro (que sabemos es la tercera parte de área de la base por la altura), de modo que el flujo puede 

calcularse efectuando la siguiente cuenta: jj s f •ñdS = 5 = - (^) 4 = ^p-, igual que antes. 



Ilustración 4. Las líneas de campo y el punto P G = (1, 1, 1) y el plano S 


Ejercicio 4. Sea C la curva frontera, orientada positivamente, del recinto D, descripto en coordenadas polares (r, 0) por las 
inecuaciones 0 < r < 2 (cos(0) + sen(0)). Definir una función escalar h: M — » M tal que la circulación de f sobre C, coincida 
con el área de D, siendo f(x, y) = Vh(x 2 y) + (xe x + 2xy,x h(x) - ey eos (y 3 ))- 


La curva C es un circunferencia centrada en P 0 = (1, 1) que pasa por el origen de coordenadas, pues r = 2 (cos(0) + sen(0)) 
sii 0 r 2 = 2 r (cos(0) + sen(0)) sii x 2 + y 2 = 2x + 2y, esto es (x-1) 2 + (y-1) 2 = 2. De los dos sumandos de f, sólo nos 
ocupamos del segundo, pues la circulación de Vh(x 2 y) es claramente nula, para cualquier h diferenciable x . De esta manera, 
la circulación resulta de la contribución del segundo miembro. Estando la curva C positivamente orientada xi , y suponiendo 
que pudiéramos hallar una h de clase C 1 en M, es claro que el campo vectorial f también será de clase C 1 en W, de modo 
que podrá aplicarse el teorema de Green y escribir entonces, llamando P(x, y) = x e x + 2 x y, Q(x, y) = x h(x) - ey cos(y 3 ) 




[' Qx (.X, y) - P y (x, y)] dx dy 



[xh\x) + h(x) — 2x] dx dy 


Es suficiente™ que h sea una función derivable en M que satisfaga la ecuación diferencial xh'(x) + h(x) - 2x =1 para 
cumplir con el enunciado. La ecuación diferencial es [x h(x)]' = 1 + 2 x, una de cuyas soluciones (la única definida en M) es 
h: M — > M tal que h(x) = x + 1. Lo verificamos: VxeM: [x h(x)]' = [x (x +1)]' = [x 2 + x]' = 2x + 1. Podemos entonces 
responder que una función que satisface el problema™ 1 es h: M — > M. tal que h(x) = x + 1. 


a 


Nuevamente, decimos que h es la única solución de la ecuación diferencial que también satisface el enunciado, pero no 
es la única función que satisface el enunciado, desde que la ecuación diferencial resultó una manera (y no la única) de 
asegurarse que la circulación y el área coincidan. Para mostrar que hay otra podríamos buscar una función hi: M — > M. tal 


x Se da por sabido que la circulación de un campo de gradientes sobre un lazo se anula, cuya prueba es directa conociendo el resultado de la independencia del 
modo en que se recorre el camino entre dos puntos. 

X1 De su regularidad no dudamos, basta la parametrización simple clásica y(t) = (1+ a cos(t), 1 + a sen(t)), con a = V 2, te [0, 2 ti), que la provee de la orientación 
necesaria, con y’ nunca nulo pues para todo t es || y’|| = a ^ 0. La Inyectividad es también directa (probarla). 

xu Y de ninguna manera necesario : es falso que debamos exigir esa condición para lograr que la circulación coincida con el área de la región. Existen infinitas 
funciones h que no cumplen la condición y sin embargo satisfacen el problema. La diferencia es que la función h que hallamos aquí, logra hacer que la circulación 
sobre cualquier lazo coincida con el área de su interior. 

xm Observamos que habíamos supuesto a priori, para aplicar el teorema de Green, la suavidad de h de primer orden. A posteriori podemos decir que la hipótesis se 
cumple con nuestra función h, dado que es C°°(M). 
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que hi (x) = k, con k a determinar de modo que ajuste el valor del área (siendo 31 un círculo de radio V2, es área (31) = 2n), 
esto es que debe cumplirse la siguiente condición: 



[xlí (x) + h(x) — 2x\ 


dx dy = 



[k — 2x\ dx dy = 2n 


El cálculo de la integral se simplifica si se observa que k-2x = k- 2- 2 (x -1), y dado que la integral de x - 1 en 31 es nula 
(es el momento estático respecto a un eje baricéntrico), sólo queda la integral de la constante k - 2, que naturalmente es el 
producto de k - 2 por el área de 31, resulatando entonces (k - 2) 2 tt = 2n, de donde k = 3. Tenemos así otra función que 
satisface lo pedido: hi: M — » M tal que hi (x) = 3. De que hay infinitas es (casi) evidente xiv . 


Ejercicio 5. Hallar los puntos de la superficie S = {(x ; y, z)et 3 : z = Jxy + 1 } más próximos al origen de coordenadas. 


El dominio de la función z: D z aW — >E es D z ={(x, y)el 2 : x y + 1 > 0}. Debemos investigar los extremos de la función 
distancia d(X, 0) = ||X||. Dado que los extremos de ||X|| y los de | |X| | 2 son los mismos xv , trabajamos con la función de tres 
variables (ligadas) f: IR 3 — »E tal que f(X) = | |X| | 2 = x 2 + y 2 + z 2 , sujeto a la restricción XeS. Pero esto es lo mismo que que 
examinar los extremos del campo escalar de dos variables (libres) g : D z aW — »E tal que g(x, y) = x 2 + y 2 + xy + 1. Para 
resolver este problema estudiamos primero los extremos en la frontera de D z , esto es en el conjunto de nivel 0 de la función 
z, C(z, 0) = dD z = {(x, y)el 2 : x y = - 1}. 

Es claro que dD z , ={(t, -1/t), t ^0, te IR}, de modo que los extremos de g en dD z , son los extremos de la función de una 
variable h(t) = t 2 + 1/t 2 , cuyo gráfico se muestra en la ilustración 4. El mínimo local de esta función escalar se alcanza en 
los puntos cuyas abscisas satisfacen h'(t) = 2t-2/t 3 = 0, esto es 1 1 = -1 o t 2 = 1 y su valor es 2. De esta manera, hemos 
encontrado que el campo escalar g alcanza el mínimo 2 en los puntos Ai = (1, -1), A 2 = (-1, 1), lo que indica que, 
regresando al problema original, entre los puntos de S en los que z = 0, los que se hallan más próximos al origen son Pi = 
(1, -1, 0) y P 2 = (-1, 1, 0), siendo tal distancia V 2. 



Ilustración 5. El gráfico de h: E -{0}— >E tal que h(t) = t 2 + t 2 

Ahora es necesario averiguar lo que sucede en el resto de los puntos de S, esto es buscamos los mínimos de g en el interior 
de D z , el conjunto abierto D = {(x, y)eIR 2 : x y + 1 > 0}, y resolviendo Vg(x, y) = (2x + y, 2y + x) = (0, 0) sólo se tiene el 

punto P 0 = (0, 0) con un Hessiano H g (Po) =det > 0' de m °do que en el punto Mo = (0, 0, 1) e S se tiene un mínimo 

de la distancia xvi . 


X1V Como ejercicio, comprobar que la función I 12 : M — > M. tal que hi (x) = 3x/2 también satisface el problema ¿podría dar una expresión explícita que muestre que, 
efectivamente, existen infinitas funciones que satisfacen el problema? 
xv ¿por qué? 

XV1 Aquí suponemos conocido el resultado de que un campo escalar g de dos variable de clase C 2 en un punto interior P G en el que se anula su gradiente tiene un 
extremo (de hecho, estricto) en ese punto si su hessiano en ese punto H(P 0 ) es positivo; el extremo es mínimo si Dng (P 0 ) es negativa. Ver Rey Pastor, J., Pi Calleja, 
P., & Trejo, C. (1968). Análisis Matemático II. Cálculo infinitesimal de varias variables. Aplicaciones (Séptima edición ed., Vol. II). Buenos Aires: Kapelusz, 
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Finalmente, es preciso decidir cuál de estos tres puntos hallados (Pi, P2 y Mo) es el que se halla más próximo. Dado que 
d(Mo, 0) = 1 < d(Pi, 0) = d(P2, 0) = V 2„ finalmente concluimos que M 0 es el (único) punto má próximo al origen de 
coordenadas, y su distancia a él es 1 . Eso es lo que ilustra el siguiente gráfico, donde se aprecian sobre la superficie los dos 
puntos Pi y P2 (color blanco) y el punto Mo más próximo al origen (color verde), que le da el mínimo a la función f. 




Tal vez sirva mencionar algo respecto de la dimensión de los espacios sobre los que se trabaja, pues parece haber 
resultado alguna confusión, a juzgar por algunas resoluciones. Es verdad que el ejercicio puede resolverse de un modo muy 
sencillo, observando lo que sucede con los conjuntos de nivel de la función a extremar y observando la relación que 
guardan con la restricción, ambos en el mismo espacio, natur^mente. El gráfico siguiente muestra las esferas (conjuntos de 
nivel, que van aumentando de tamaño hasta contactar por pernera vez a la superficie restricción S en el punto M 0 : allí está 
el mínimo, y asunto terminado. 



§70.2, Teorema 1, p.205, o también Apóstol, T. (1980). Calculus volumen 2. Cálculo con funciones de varias variables y álgebra lineal, con aplicaciones a las 
ecuaciones diferenciales y a las probabilidades (Segunda edición en castellano [Original: Calculus II, Multi -variable calculus and linear algebra, with applications 
to differential equations and probabilty] ed., Vol. II). (F. Vélez Cantarell, Trad.) Barcelona: Reverté, §9.12, Teorema 9.7, p.380 

El grado de regularidad suficiente para la validez de la afirmación es C 2 , sin embargo, en algunos textos se pide C 3 , lo que facilita la prueba; ver por ejemplo, 
Curtis, P. (1979). Cálculo de varias variables con álgebra lineal (Primera edición. Primera reimpresión [Original 1970: Multivariate Calculus with Linear Algebra] 
ed.). (M. C. Sangines de Salinas, & O. Mourut de Montpellier (revisora), Trads.) México D. F.: Limusa., Teorema 7.7, p. 329; Marsden, J. E., & Tromba, A. J. 
(1991). Cálculo Vectorial (Cuarta edición en español del original Vector Calculus, Third Edition. ed.). (M. López Mateos, & S. Adarve, Trads.) Wilmington, 
Delaware, EUA.: Addison-Wesley Iberoamericana, §4.2, Teorema 4, p.253. El teorema para un campo de dos variables, con regularidad C , se encuentra en un 
artículo de Lagrange de 1759, titulado Recherches sur la méthode de maximis et minimis, Cf. Hairer, E., & Wanner, G. (1996). Analysis by its History (Primera 
edición ed.). New York: Spinger-Verlag.Chapter I, Higher Derivatives and Taylor Series, Theorem 4.7, p. 324. La generalización a varias variables introduce el 
actual hessiano en el artículo de Otto Hesse de 1857: Ueber die Criterien des Máximums und Minimus der einfachen Intégrale. 
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En este documento encuentran referencias de los siguientes textos, de los que sería deseable que cada uno de los alumnos tuviese, 
ciertamente, al menos uno. En las clases hemos hecho continua referencia al texto de Apóstol y al de Marsden-Tromba. 

CQ Apóstol, T. (1980). Calculus volumen 2. Cálculo con funciones de varias variables y álgebra lineal, con aplicaciones a las 
ecuaciones diferenciales y a las probabilidades (Segunda edición en castellano [Original: Calculus II, Multi -variable calculus and 
linear algebra, with applications to differential equations and probabilty] ed., Vol. II). (F. Vélez Cantarell, Trad.) Barcelona: 
Reverté. 

CQ Curtís, P. (1979). Cálculo de varias variables con álgebra lineal (Primera edición. Primera reimpresión [Original 1970: 
Multivariate Calculus with Linear Algebra] ed.). (M. C. Sangines de Salinas, & O. Mourut de Montppellier (revisora), Trads.) 
México D. F.: Limusa. 

CQ Hairer, E., & Wanner, G. (1996). Analysis by its History (Primera edición ed.). New York: Spinger-Verlag. 

CQ Lang, S. (1990). Introducción al Análisis Matemático (Primera edición. [Original 1968, Undergrauate Analysis] ed.). (M. 
López Mateos, & M. Muñoz Mella, Trads.) Wimington, Delaware.: Addison Wesley Iberoamericana. 

CQ Marsden, J. E., & Tromba, A. J. (1991). Cálculo Vectorial (Cuarta edición en español del original Vector Calculus, Third 
Edition. ed.). (M. López Mateos, & S. Adarve, Trads.) Wilmington, Delaware, EUA.: Addison -Wesley Iberoamericana. 

CQ Rey Pastor, J., Pi Calleja, P., & Trejo, C. (1968). Análisis Matemático II. Cálculo infinitesimal de varias variables. 
Aplicaciones (Séptima edición ed., Vol. II). Buenos Aires: Kapelusz. 

CQ Santaló, L. (1993). Vectores y tensores con sus aplicaciones (Décimocuarta edición ed.). Buenos Aires: Eudeba. 
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